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Abstract. We apply EROS photometric data to interpret previously published Keck 
and VLT spectra of the binary-microlens caustic-crossing event EROS-BLG-2000-5. We 
show that the VLT data imply that the outer ~ 4% of the limb of the K3-giant source is 
strongly in emission in Ha, in contradiction to available models of the photosphere. This 
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conflict could be resolved if the integrated Ha emission from the chromosphere were equal 
to 8% of the integrated Ha absorption from the source as a whole. These conclusions 
regarding the extreme limb are almost completely model-independent. We also present 
a general method for using the photometric data to facilitate direct comparison between 
the atmospheric model and the spectroscopic data. While this method has some model- 
dependent features, it is fairly robust and can serve to guide the analysis of spectra while 
more detailed models of the lens geometry are being developed. In addition, we find that 
the color of the limb of the source (outer 5.5% by radius) is A{V — /) ~ 0.37 redder than 
the source as a whole, so that it has the color of a MO giant. 



Key words: gravitational lensing - techniques: high angular resolution techniques: spec- 
troscopic - stars: atmospheres 

1. Introduction 

Microlensing provides a potentially powerful probe of stellar atmospheres. If a lens caustic 
(region of formally infinite magnification) passes over the face of the source, then different 
parts of the atmosphere become highly magnified at different times. By combining a time 
series of photometric or spectroscopic observations, one can therefore hope to deconvolve 
the spatial structure of the atmosphere. For point lenses, the caustics are point-like, and 
therefore the probability of such a caustic crossing is small. However, if the lens is a 
binary with components separated by of order the Einstein radius, the caustics form one 
to three concave polygons whose total length is of order the Einstein radius, and therefore 
the probability of a caustic crossing is much larger. 

Intensive photometric observations of four such binary caustic crossing events have 
yielded limb- darkening measurements (Afonso et al. 2000; Albrow et al. 1999,2000,2001a). 
A more ambitious project would be to obtain a similar intensive series of spectroscopic 
measurements. Since light from the limb of the star originates higher in the atmosphere 
than light from the center, deconvolution of a set of spectral measurements effectively 
resolves the atmosphere as a function of height (Valls-Gabaud 1998; Heyrovsky, Sasselov 
& Loeb 2001; Gaudi & Gould 1999). Until the advent of microlensing, such spatially 
resolved spectra had been obtained only for the Sun. 

Early spectra of caustic-crossing events were obtained by Lennon et al. (1996) and 
Alcock et al. (1997), but no clear differences were detected between spectra taken at 
different times. 
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The binary caustic crossing event EROS-BLG-2000-5 has now yielded the first big 
breakthrough. This event was alerted by EROSQ on 5 May 2000. On 8 June 2000, MPS0 
issued an anomaly alert, saying that a caustic crossing was in progress. Subsequent in- 
tensive observations by PLANET^ allowed them to predict the time of the second caustic 
crossing and, very importantly, that this crossing would last an unusually long 4 days. 
This prediction permitted two groups to acquire spectra with large telescopes on several 
successive nights: Castro et al. (2001) obtained HIRES {R ^ 40,000) spectra using Keck 
on the last two nights of the crossing, while Albrow et al. (2001b) obtained low-resolution 
(i? ^ 1000) FORSl spectra using VLT on all four nights. In addition, Albrow et al. 
(2001b) obtained a spectrum well before the caustic crossing, when the magnification of 
the source was approximately uniform. 

To date, each group has published only the Ha lines, but both sets of spectra cover 
several thousand A, and could potentially yield a wealth of information about the K3 
giant source. 

Here we use EROS photometric observations of this event to aid in the interpretation 
of the spectral data that have been published. In § 2, we present the EROS photometric 
data of the event. In § 3, we present an essentially model-independent argument to show 
that the VLT data imply that the outer ~ 4% of the limb of the source is very strongly 
in emission in Ha. In § 4, we develop a general method to apply photometric data of 
a caustic crossing to facilitate comparison between spectroscopic data and atmospheric 
models and apply this method to EROS-BLG-2000-5. In § 5, we characterize the limits 
of our approach. In § 6, we measure the color of the source limb by applying a slightly 
modified version of our approach to the EROS photometric data. Finally, in § 7, we 
speculate on the possibility that the Ha emission detected from the limb of the source is 
actually due to the chromosphere, and quantify the strength of chromospheric emission 
required to explain the observations. 

If the geometry of EROS-BLG-2000-5 were perfectly modelled, then there would 
be no need for the techniques introduced in this paper. However, the event is quite 
complex and difficult to model, so additional techniques are required to begin quantitative 
investigation of the spectroscopic data while a better geometric model is being developed. 
Moreover, it is likely that in future events, one will face a qualitatively similar situation, 
and the methods presented in this paper will be applicable long before the event as a 
whole is analyzed at a satisfactory level. 



^ http:/ /www-dapnia.cea.fr/Spp/Experiences/EROS/alertes.html 
^ http:/ /bustard. phys.nd.edu/MPS/index. html 
^ http://thales.astro.rug.nl/ planet/ 
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2. Photometric Data 

EROS observations were carried out using the 1 m Marly telescope at La Silla, Chile, in 
two bands Ve and Ie- EROS Ve is centered midway between Johnson V and Cousins 
R, while EROS Ie is similar to Cousins /, but broader. The transformation between the 
Eros system (Ve, Ie) and the standard Johnson-Cousins {V,I) system was determined 
by observing Landolt (1992) standards near the galactic poles and Paczynski et al (1999) 
secondary standards in the Baade window: 

Ve = 0.71 y + 0.29/ + const, /e = -0.02 V + 1.02 /+ const, (1) 

These equations taken together together imply A(Ve — Ie) — 0.73A(y — /). The first 
of equations ([T]) implies that the centroid of the Ve band is approximately 6300 A. We 
therefore present mainly Ve data because its centroid is close to Ha. EROS observations 
are normally carried out in survey mode, but owing to the importance of this event we 
devoted all available time to it during the caustic crossing. Generally, the weather was 
good, but on the final night there were intermittent clouds, and then finally we were 
clouded out entirely four hours before dawn. 

The data were initially reduced on site using our standard PEIDA (Ansari, et al. 
1996) PSF-based software, but were then reprocessed using the ISIS (Alard, 2000) image 
subtraction program. We then determined the zero point of the image-subtracted pho- 
tometry by fitting for the mean offset between it and the PEIDA photometry. Figure |l] 
shows the naive magnification Ay , which we determine by dividing the measured flux by 
the baseline flux. The curve will be explained in § 4, but from the data points alone, it is 
clear that the caustic crossing ended some time before JD' = 1733.6506: the light curve is 
almost perfectly fiat commencing at this time and continuing for 58 minutes, with a mea- 
sured slope (combining information from both bands) of d\nA/dt ~ 0.11 ± 0.08 day~^. 
(Here JD' ==JD-2450000.) 

The crosses indicate the times of the Keck and VLT observations. Note that the 
VLT observations were all coincident with our photometric observations, while the Keck 
observations both took place about 2.2 hrs after the end of our night. 

3. Ha in Emission? 

The last VLT equivalent width (EW) measurement, which was taken just before the 
source finally exited the caustic, is much lower than any of the other EW measurements. 
As we will show a few paragraphs below, this can be interpreted as indicating strong Ha 
emission from the edge of the star. In order to quantify this conclusion, we must first 
derive some relations between the photometric and spectroscopic data. 
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Let ^(u) be the magnification field of the binary as a function of angular position 
u normahzed to the Einstein radius, let £'(u) be the equivalent width of the Ha line, 
and let S{u.) be the source surface brightness in the neighborhood of this line. We will 
split the magnification field into two components A{u) = A2{u) + A3(u) corresponding 
respectively to the two images whose magnification diverges (and then vanishes) at the 
caustic, and the three non-divergent images. The observed magnification is then, 
, f (PuA(u)S(u) ^ , , , 

where 

A2(u,) = I Y^^' J ^ 

and where we have assumed that the non-divergent magnification field can be replaced by 
its value at the center of the source Uc. This is strictly true whenever this magnification 
field is either flat or has a uniform slope. We will discuss in § 5 why we believe this to 
be an extremely good approximation in the present case. 
Similarly, the observed Ha EW will be 

_ J d^uAiu)S{u)Eiu) _ A3iu,)Eo +M^c)M^c) , . 

^"^^ Jd^uAin)Siu) A3K) + A^K) ' ^' 

where, Eq is the EW of the unmagnified source, and 

_ J d^uA2{u)S{u)Eiu) 

J^2(Uc) = -=- ^ (5) 

is the EW that would be seen if a separate spectrum could be taken of the two highly 
magnified images. Thus, if Eq and E{uc) are measured spectroscopically and if the frac- 
tion of light coming from the highly magnified images /(uc) = A2{uc) / A{uc) is known 
from the photometric light curve, then one can invert equation (^) to determine £'2(uc), 

^ = 1 + ^, A.^, (6) 
Eo f Eo ' ■' A ^ ' 

If we assume no blending (see § 5), then A{uc) is directly measured from the light 
curve (Fig. |l|). In principle, A3 is changing with time and so can only be known if there 
is a full model of the light curve. However, in general A3 is expected to change slowly. 
Moreover, in the present case it is measured beginning at JD' = 1733.65, only 2.5 hours 
after the last VLT measurement (JD' = 1733.54) and is found to be nearly constant for 
an hour after that. We will therefore assume that A3 is constant throughout the crossing 
and will discuss possible small corrections to this approximation in § 5. 

For the last VLT point, Albrow et al. (2001b) report !E(1733.54) ^ 0.79 ± 0.03 A, 
and they adopt Eq = 1.03 A, i.e. A = -0.23 ± 0.03. From the EROS lightcurve (Fig. |l|), 
/(1733.54) = 0.060 ± 0.007. As we will show in § 4, at this time the source center was 
approximately rj ~ 0.963 source radii past the caustic, so that i?2 (1733.54) is a weighted 
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average of the outer 4% of the source. Here rj is the perpendicular distance from the 
source center to the caustic in units of the source radius, and is taken to be negative 
when the source center is inside the caustic. Combining these facts with equation (^) 
yields 

^(11?^ = -2.9 ±0.7 (,.0.96), (7) 

-C'O 

which is to say that the limb is in emission more strongly than the star as a whole is in 
absorption. As we show in § 4, this would strongly contradict the model of Albrow et al. 
(2001b). 

We note that since the caustic is a square-root singularity (and assuming circular 
symmetry of the source) , 

drrS(r)E(r) 

^-Vt^ (8) 
drrS{r) 

where r is normalized to the source radius. That is, E2 is very nearly equivalent to a 
simple weighted average of the EW by the light coming from the entire periphery of the 
source. 

Applying equation (ph to all the VLT and Keck measurements, we find. 



Table 1. Implications of Measured Ha EWs 



Telescope 


JD' 


E{A) 


aiE) (A) 


/ 


^(/) 


V 


£^ (A) 


aiE,) (A) 


VLT 


1730.60 


1.06 


0.01 


0.846 


0.001 


-0.595 


+ 1.07 


0.01 


VLT 


1731.67 


1.09 


0.02 


0.803 


0.001 


+0.097 


+ 1.10 


0.02 


VLT 


1732.66 


0.98 


0.01 


0.626 


0.001 


+0.611 


+0.95 


0.02 


VLT 


1733.54 


0.79 


0.03 


0.060 


0.007 


+0.963 


-2.97 


0.68 


Keck 


1731.953 


0.944 


0.002 


0.773 


0.001 


+0.256 






Keck 


1732.950 


0.869 


0.006 


0.510 


0.002 


+0.738 







Note that since Castro et al. (2001) did not measure Eq, we cannot use equation 
(^ to evaluate E2, for the Keck data. One cannot simply adopt the VLT value of Eq 
because the VLT and Keck data are affected by systematic differences as discussed in the 
penultimate paragraph of § 4. We note that D. Minniti (2001, private communication) 
has high-resolution VLT data for this event when it was highly (but not differentially) 
magnified, so it may eventually be possible to correct this shortcoming. 

The results reported in Table 1 can be directly compared with results from model 
atmospheres, subject to some minor qualifications discussed in § 5. 
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4. Predicting Equivalent Widths from the Lightcurve 

Without loss of generality, one can write the magnification as 

^(u) = A3(u) + Z(u)Au;'/'e(Aux), (9) 

where Au± is the perpendicular distance from u to the caustic and is a step function. 
As we discussed earlier, unless the source is actually probing the cusp region, it is an 
excellent approximation to substitute A3(u) ^3(uc). Since the cusp approach occurred 
^ 4 days after the caustic exit, and since ^3(uc) was observed to be almost perfectly flat 
for an hour after the crossing, we adopt this approximation. For closely related reasons 
(discussed in § 5), we adopt Z(u) Z{uc). Under these assumptions, and making use 
of the fact (explicitly demonstrated in § 5) that the caustic curvature is small on scales 
of the source, the spatial averaging over the model atmosphere (required to compare 
theoretical predictions with the empirical results summarized in Table 1) is independent 
of Z, 



^2(7?) ^ Z(u,)G(r;), G(7?) = ^ /' dx{x - ry)-!/^ dyS{r), (10) 



E2{ri) - i?(r/) - ^ f dx{x-ir"^ ,dyE{r)S{r) (11) 



G(r7) ' 



where x and y are dimensionless variables and r = \/ x"^ + J/^- Note that for a uniform 
source, G{rj) is given by e.g., Figure 1 of Gould & Andronov (1999). 

In this section we will make the additional assumption that ^3 and Z are independent 

of Uc, 

^3(uc) Ai, Z{uc) — > Z (assumption). (12) 

As we discuss in § 5, this assumption is less robust than others made in this paper, but 
is still quite reasonable for our purposes. 

Albrow et al. (2001b) have presented calculations of A{'ri) and E{r]) for a K3 giant 
atmosphere model in their Figures 1 and 4 respectively, under the assumption that equa- 
tion (|l2|) is valid, for two specific choices of the pair (^43, Z). The resulting curves were 
intended only for schematic purposes, and in particular assumed that 77 was a linear 
function of time [7/ = (JD' — 1731.8)/1.9], which is known not to hold in the present case 
(see eq. |l^, below). However, our interest in these curves is solely that they allow us to 
evaluate the functions G and H that derive from the underlying atmosphere model. As 
a matter of practical computation, we must normalize G, and so we adopt Gmax ~ 1.40. 
However, the final predictions of EW are completely invariant under changes of this 
parameter. 

We then apply these extracted functions to the EROS data. First, we measure A3 = 
3.2294 ± 0.0034, where the error includes only the measurement error of the post-caustic 
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flux (and not the baseline flux that normalizes it). This is because all predictions are 
invariant under changes of the assumed baseline flux. We also measure Amax ~ 21.02 
from the peak of the curve, and hence find Z = (A,„ax — ^3)/Gmax — 12.7. Then for 
each observed magnification, A, we determine — G^^[{A— A^j/Z]. Since G^^ is double 
valued we must assign a branch depending on whether JD' is before or after 1730.59 
(the time of the observed peak of the light curve), and for points that have A > Amax, 
due to observational error, we assign ?/ = G^^[(A,nax — ^s)/^] - The resulting values of 77 
are shown in Figure For most caustic crossings, ri{t) is reasonably approximated by a 
straight line, but clearly this is not so in the present case. Indeed from the cusp approach 
(sharp peak) in the light curve four days after the caustic crossing, it is known that the 
source must have exited nearly tangent to the caustic, so that the caustic curvature must 
play an important role. To interpolate to times when we have no data, we fit therefore 
the points 77 (i) to a parabola whose equation is given by, 

77 = 1 + 0.3298(JD' - ieccnd) - 0.0637(JD' - tcc,ond)', tccond = 1733.6506, (13) 

and whose form is shown in Figure |^. This curve is a good fit to the data except for the 
first night where it deviates by approximately the scatter in the points. This shows that 
at least one of our assumptions is not fully satisfied: either the caustic deviates from a 
parabolic form, or the assumption (^2|) is too much of an oversimplification, or perhaps 
the limb-darkening model of Albrow et al. (2001b) does not perfectly reproduce G{ri) 
for the EROS Ve band. Whatever the cause of the deviation, the fact that it is more 
pronounced on the first night is simply due to the fact that the errors in 77 are larger. 
This in turn can be traced to the fact that this night covers the peak of the curve (see 
Fig. I), where dG/drj - 0. The error in the estimate of 77 for each measured point is given 
by a{ri) = (T{Fv)/{Fv.hasc\incZ\dG/dri\). These larger errors are reflected in the larger 
scatter of the points on the first night in Figure ^. The errors are not displayed in Figure 
^ in order to avoid clutter, but are shown in a different projection in § 5. 

Note that the error on the time of the end of the caustic crossing (77 = 1), derived from 
the fit leading to equation (|l3|), is very small, tec, end — 1733.6506 ± 0.0048. However, the 
actual error is dominated by the assumption that the limb-darkening model of the curve 
shown in Figure |l| (ultimately derived from Fig. 1 of Albrow et al. 2001b) is correct, which 
it almost certainly is not at this level. Nevertheless, ice, end cannot be much after this best- 
fit value because the observed light curve is flat (or very slightly rising) beginning at this 
time. 

Finally, we predict the EW, E{t), using the values of ri{t) found from equation (p^), the 
tabulated functions G{ri) and H{rj), the measured values of A(t) and A3, and equations 
(^ and (|l^) . The result is shown in Figure ||. 



EROS Collaboration: Contraints on EROS-BLG-2000-5 Spectroscopy 9 

As can be seen from Table 1, the VLT EWs are about 0.1 A larger than the corre- 
sponding Keck EWs, even taking account of the fact that the observations were a few 
hours earlier. One possible explanation for this discrepancy is that the low-resolution 
VLT spectra were affected by blended lines that were removed from the Keck HIRES 
spectra. From Figure 2 of Castro et al. (2000), this indeed appears plausible. For purposes 
of displaying our results in Figure ^, we therefore assume that 0.1 A of the VLT EWs are 
due to blends. Hence, we reduce all of them by this amount and correspondingly rescale 
the Albrow et al. (2001b) model downward by 10%. There remains the freedom to adjust 
the relative offsets between the Keck and VLT points (because our 10% VLT blending 
estimate is only a guess), and to adjust the overall scale of the curve (because this is set 
by Eq which is measured from the VLT pre-caustic spectrum - but with a large error 
bar). 

Nevertheless, there is no amount of adjustment that can reconcile the points with the 
curve. The point with the largest deviation between observed and predicted EW is the 
last VLT spectrum. In the previous section we showed that the model could predict such 
a low EW only if the outer 4% of the source were strongly in emission. By contrast, the 
model of Albrow et al. (2001b) predicts that the outer 4% will have substantial (albeit 
reduced) absorption: £^(0.963) = i7(0.963)/G(0.963) = 0.0105 A/0. 0184=0.57 A (with 
no correction for line blends). 

5. Examination of Assumptions 

We have used EROS lightcurve data to facilitate the comparison of spectra taken during 
the caustic crossing of EROS-BLG-2000-5 with an atmosphere model of the K3 source. 
An alternative approach would be to fit the lightcurve data to a lens model and then 
apply the lens model directly to the atmosphere model to make predictions about the 
spectra. The latter approach would appear to be more secure but, unfortunately, to date 
there are no published models of this event, and the event may prove difficult to model. 
Thus, since our approach is the only one available at present, it will be worthwhile to go 
through the assumptions that were introduced and see what order of errors they induce. 
This review can also serve as a checklist if our method is applied to future events. 

We assumed no blending. There are three lines of evidence that this is a good assump- 
tion. First, the lensed source is a bright giant, so it is fairly improbable that it would 
be seriously blended. Second, there is no evidence for blended light in the spectrum of 
Castro et al. (2001) as there would be if there were a bright blend offset by at least 
15kms~^ in radial velocity from the source. Third, the event is nearly achromatic (see 
Fig. This means that either the blend is very weak in which case it has little effect, 
or it has virtually the same color as the source, and so very likely has the same Ha EW 
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as well. To the extent that this is the case, the predicted EWs are exactly the same as 
in the unblended case. 

Next, we assumed that the non-divergent magnification field ^3(u) could be replaced 
by its value at the center of the source, ^3(uc). This is valid to the extent that on 
scales of the source size, the field is equal to a constant plus a gradient. The only 
way this could fail to hold would be if the source were probing the neighborhood of a 
cusp. It is clear that the source is approaching a cusp, because the lightcurve betrays a 
characteristic cusp-approach behavior 4 days after the caustic exit. However, in the hour 
after the caustic exit, the light curve is nearly fiat, indicating that this cusp is not yet 
significanly affecting the magnification pattern at the time of the exit. 

Third, we assumed that A3 would be constant throughout the crossing. The measured 
logarithmic slope of A3 immediately after the crossing is dlnA^/dt = (0.11±0.08) day^^, 
so that at the time of the last spectrum measurement (0.11 days before the end of the 
crossing), A3 was within a percent or two of the final value. A 1% change in A3 leads 
to a change of ^ —17% in / and so of ~ —0.5 in equation (^. This is not negligible, 
but it would not qualitatively affect the conclusion that the outer 4% of the source is 
in emission. Plausibly, at much earlier times ^3 could have differed from its final value 
by ~ 10%, but as these times such a change would have a very small impact on the 
conclusions simply because A3 does not dominate the total fiux. For example, suppose 
this had been the case at JD' 1732.66. Then / would have fallen from 0.626 to 0.588, 
which in turn would change E2 by only 0.005 A. These are all the assumptions needed 
through § 3. 

In § 4, we went on to make three additional assumptions. The first was that on 
the scale of the source the caustic could be approximated as a straight line. From the 
measurements of the first crossing, the source is known to cross its own diameter in 
^ 1 day. If we interpret the quadratic behavior of equation ( p!3| ) as arising from the 
curvature of the caustic, then over a diameter, the caustic deviates from a straight line 
by ~ 8 X 10^^ of a source radius, certainly small enough to ignore. 

Second, we assumed that at any instant Z{u) could be replaced by its value at its 
center Z{uc). Or equivalently, that we are in the square-root singularity regime and the 
coefficient of this regime varies at most linearly over a source diameter. From Table 1 
and the adopted value of Z = 12.7, one sees that that at the trailing limb of the source 
for the first measured point, A2{u) = 10, which is normally well into the square-root 
singularity. Since the source diameter is ^ 1 day and at the exit it is still 4 days from 
the post-caustic peak in the light curve due to a cusp approach, the assumption that 
quadratic variation of Z is not significant on the scale of the source size appears very 
reasonable, but by no means proven. 
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Third, we assumed that ^3 and Z remain constant throughout the crossing. We 
discussed above why this is a reasonable assumption for Aj,. The argument regarding Z 
is similar. Suppose that Z in fact fell linearly by 20% from the time of the peak until 
the end of the crossing (and suppose ^3 remained constant). The product G{rj)Z = A2 
is then an observable and so remains fixed. Thus, at the last measurement, Z would be 
19% lower, so G{ri) would actually be 24% higher than we have allowed. This would drive 
r] down from 0.963 to 0.956, but it would affect the EW prediction shown in Figure ^ by 
well under 1%. Hence, while it would be better not to have to make this approximation, 
the errors it induces are extremely small compared to the level of disagreement between 
the models and the data in Figure ^. 

To directly probe the scale of the systematic errors introduced by our approxima- 
tions, we conduct the following test. Instead of adopting the Gp{ri) derived from the 
atmosphere model of Albrow et al. (2001b), we directly fit for Gv'(?7;c) from the ob- 
served Ve magnification, Av(t), using the 6-parameter functional form, 

Av{t)^A3 + ZGv[v{t);c], (14) 

where c is a linear limb-darkening coefficient, and rj is a (3-parameter) quadratic function 
of time as in equation (p^. We find c ~ 0.644±0.032, tec. end — 1733.6280, and linear and 
quadratic coefficients 0.3421 and —0.0596, respectively. The resulting curve ri{t) (Fig. ^) 
is very similar to the corresponding curve (Fig. ^ for Gp. Moreover, when we evaluate rj 
using Gv instead of Gp, we find that the predicted EWs are almost exactly the same as 
those shown in Figure ^: the changes are an order of magnitude smaller than the error 
bars on the VLT and Keck data points. 

Figure ^ shows the photometric residuals of the fit during the first night (JD' ~ 
1730.6), for which the residuals in terms of 77 seen in Figure ^ are most pronounced. 
Typically, these residuals are extremely small, < 0.005 mag, and are of order the individ- 
ual errors. However, taken together they show a clear trend which is highly significant: 
the observed light curve peaks earlier than the model. The offset is about 0.1 day, or 
A77 ^ 1/15. This level of discrepancy induces an offset in the predicted EW on the first 
night of about 1/2 the statistical error of the EW measurement (see Fig. However, 
since the model is only used to interpolate to times when there are no photometric data, 
and since the model and data are in excellent agreement in the neighborhoods of the 
Keck points (which are the only ones where there is no coincident photometry), this 
discrepancy has no practical impact. 

Parenthetically, we note that while one might be tempted to take the good agreement 
shown in Figure |^ as supporting the believability of the derived limb-darkening coefficient, 
in fact this determination must be viewed skeptically. Any change in A3 or Z during the 
very long crossing, as well as any deviation of the caustic from a purely parabolic form. 
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will be subsumed in the limb-darkening coefRcient. Hence, a valid measurement of c 
can be obtained only in conjunction with a complete determination of the geometry 
of the binary lens. The main implication of this fit is that, while the limb-darkening 
coefficient may be significantly affected by these various deviations from a simple model, 
the predictions of equivalent widths are not. 

We conclude that all the assumptions we have made are quite reasonable, at least for 
the present case. Future applications of this method will have to be carefully justified 
based on an analysis of the event that is parallel to the one given here. 

6. Color of the Limb 

Essentially the same method presented above can be used to find the color of the two 
highly magnified (caustic) images. Recall from equation (||) that for 1 — 77 <C 1, this is 
very nearly equivalent to the color of the annulus of the source 1 — 77 < r < 1. The 
difference in the color of the two highly magnified images relative to the unmagnified 
source is given by, 

A{Ve ~ Ie) = {Ve - Ie)2 - {Ve - /B)source = 2.5 logf ^) (15) 

\Fv,2 dFi J 

where F/.2 and Fv,2 are the fluxes due to the two caustic images, and dFy/dFj is the 
slope shown (normalized to the baselines) in Figure ^. Since this figure is restricted to 
points away from the caustics and cusps where the source is not differentially magnified, 
— 2.51og((iFv/(iF/) is the color of the source. Under the assumption (p^), 

Fv,2 = Fy — Fv.oYvd, (16) 

where Fy is the observed Ve flux at a given time and Fy^cnd is the flux just after the 
end of the caustic crossing (and similarly for Ie)- Note that equations ([l5| ) and (|l^) are 
valid independent of whether the source is blended, and that when combined they can be 
evaluated from directly observed quantities without any modeling. The only assumption 
is that A3 does not change significantly between the time of the measurement and the 
end of the caustic crossing. 

Figure ^ shows the result of applying equations ( p^ ) and (^6|) to the EROS data. 
The time coordinate has been changed to source position relative to the caustic (77) 
using equation (^3|). The color difference is small to negligible on the first two nights 

[A(Ve-/e) = 0.0006±0.0006 and -0.0015±0.0006 for 77 0.6 and 77 - 0.1 respectively] 

reflecting the fact that the caustic images receive significant contributions from a range 
of radii at these positions (Gaudi & Gould 1999). On the third night, the caustic images 
are significantly redder, A(Ve — Ie) — 0.046 ± 0.001 at 77 ^ 0.6. Here, the analog of 
equation (||) approximately holds, so that the color of the caustic images is the same 
as that of the outer 1 — 77 ^ 40% of the source (by radius). Finally, on the last night. 
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we find A{Ve — Ie) — 0.27 ± 0.02 at ry ~ 0.945, equivalent to an annulus comprising 
the outer 5.5% of tlie source. This is roughly equivalent to /S.{V ~ I) ^ 0.37 in standard 
Johnson/Cousins bands. Hence, assuming that the K3 source has a typical (y— /)o ^ 1.36 
color (Bessell & Brett 1988), the Hmb has {V - I)o ^ 1.73, i.e. the color of a MO giant. 



7. Probing the Chromosphere? 

The color (and thus temperature) difference of the source periphery relative to the source 
as a whole (A(y ~ I) ^ 0.37) is modest compared to the dramatic shift from Ha absorp- 
tion to emission implied by equation That is, MO stars absorb rather than emit Ha. 
Hence, to explain this shift it would appear necessary to invoke some additional structure 
in the source atmosphere. The chromosphere presents itself as an interesting candidate. 
We therefore evaluate what strength of chromospheric emission is required to account 
for the sharp decline in EW reported by Albrow et al. (2001b) on the last night. 

Seen in projection above the limb of the photosphere, the solar chromosphere emits 
strongly in Ha and is essentially black in the neighboring continuum. If we assume that 
the same is true of the source of EROS-BLG-2000-5, then equation (|) can be generalized 
to 



— ^3(Uc)£^0 + A2(Uc)£;2(Uc) - A2,chrom(Uc)£'o^ „^ 

where Eq still represents the EW of the unmagnified source, including the effect of Hq 
chromospheric emission, X is the ratio of this chromospheric emission to the algebraic sum 
of the Hq absorption by the stellar atmosphere and of the chromospheric emission, and 



^2,chrom(uc) is thc magnification of the chromosphere by the two caustic images. Note 



that A2,chrom IS the ratio of the H^ chromosphere photon flux in the two caustic images 
to the Hq chromosphere photon flux normally received in the absence of magnification, 
and so is not exactly analogous to A2, which is normalized to the number of photons 
normally received from the source as a whole. 

It is straightforward to show that in the limit 77 ^ 1, 

A2,chrom ^ ^ iV ^ ■ (18) 

See for example Figure 1 from Castro et al. (2001). Here, is the magnification due to 
the two caustic images of a point source lying exactly one source radius inside the caustic. 
We will adopt = 12.7, which is the value measured in § 4. As discussed in § 5, Z is 
actually a function of Uc and so in principle could be different at the peak of the light 
curve (77 ~ —0.5) where it is measured, from its value at 77 = 0.963. We therefore assign 
this parameter an error of 20%. (In addition, an error could be introduced if Gmax for 
the actual limb-darkened profile differed from the one adopted. However, for a reasonable 
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range of limb-darkening profiles, this uncertainty is < 2%.) Solving equation ( p^ for X 
and substituting in values for the various parameters, we find 



where we have taken E2/Ea — 0.57 as evaluated in § 4 based on the Albrow et al. 
(2001b) model. Note that since the ratio of the two terms within the brackets is ~ 0.05, 
uncertainties in the estimate of E2/EQ are not likely to significantly impact the final 
result. 

This value of X is about two orders of magnitude higher than the corresponding 
quantity for the Sun, but giant stars may have a more significant chromosphere than the 
Sun. Alternatively, the errors in the Ha measurement may have been underestimated. In 
any event, this estimate of X can be directly compared with chromosphere models. 
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Fig. 1. Magnification Ay of a K3 giant during the second caustic crossing of EROS- 
BLG-2000-5. Open circles show the fluxes in the EROS Ve band divided by the baseline 
flux (i.e., assuming no blending.) Large crosses represent the times of spectroscopic ob- 
servations by Albrow et al. (2001b) at the VLT (V) and by Castro et al. (2001) at Keck 
(K). The curve is the projection of the simple quadratic model shown in Fig. |^ into (mag- 
nification, time) coordinates using the stellar-profile function Giji) [eq. (p^j derived from 
Fig. 1 of Albrow et al. (2001b). (Irregularities are due to the difficulty of transcribing 
the figure.) The insert also shows EROS Ie points (crosses) for the last night. From the 
two bands together, the slope of the magnification during the hour after the caustic exit 
at JD' = 1733.65 is d\nA/dt = (0.11 ± 0.08) day^. Here, JD' =JD-2450000. 
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Fig. 2. Photometric determination of the perpendicular distance rj from the source center 
to the caustic in units of the source radius. First, we find G{ri) — {A — ^3)/^ where A 
is the observed magnification (Fig. |l|), A3 = 3.23 is the post-caustic magnification, and 
Z = {Anm^ — A3) /GniayL — 12.7 (scc tcxt). Wc then invert the stellar-profile fmiction G{ri) 
from Fig. 1 of Albrow et al. (2001b) to determine rj. The curve is a quadratic fit to all 
the points. The fit time of its endpoint (JD' = 1733.65 at 77 = 1) occurs just before the 
start of the last hour of data which, from Fig. ^, is independently known to occur after 
the caustic exit. 
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Fig. 3. Measured and predicted Ha equivalent widths (EWs). Small points show the 
predictions from the photometric data. For each measured magnification (Fig. |l|), the 
source-caustic separation 77 is determined (Fig. H), and so also the EW of the two highly 
magnified images inside the caustic, E2 = H{rj) / G{rj). Here G and H are functions 
extracted from the atmosphere model of Albrow et al. (2001b). See text. The predicted 
EW is then the weighted average of E2 and Eq, the EW of the unmagnified source. 
The curve results from applying the same calculation to the curves in Figs. |^ and ^. 
Solid points are measured EWs from VLT FORSl low-resolution data of Albrow et al. 
(2001b) and open points are from HIRES KECK data of Castro et al. (2001). Crosses 
show where these points should be if they agreed with the model. The VLT points and 
the corresponding model have been rescaled assuming a 10% blend from a neighboring 
line. This is plausible based on Fig. 2 from Castro et al. (2001) and brings the VLT points 
approximately into line with the Keck points. However, no amount of rescaling would 
bring the atmospheric model into agreement with all the data points. 
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Fig. 4. Comparison of fliixes in two bands Ve and Ie for all points on the EROS light 
curve except during the two caustic crossings and the cusp approach (when differential 
magnification could induce color changes). The regression (solid line) passes extremely 
close to the origin, showing that either the event is not significantly blended or if it is, 
the color of the blend is extremely close to that of the source. The large crosses show the 
10 outlier points (out of 250) that were excluded from fit. 
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Fig. 5. Filled points and upper curve: Photometric determination of the perpendicular 
distance rj from the source center to the caustic in units of the source radius, obtained 
by directly fitting the EROS Ve data to a limb-darkened profile function Gy (77; c), where 
the linear limb-darkening parameter c is treated as a free parameter and is found to be 
c = 0.644 ± 0.032. Open points and lower curve: Same as Fig. ^ which was based on 
the stellar profile function Gp{ri) taken from Albrow et al. (2001b). The a priori model 
Gp{ri) yields almost exactly the same curve as the limb-darkening fit GviVj c). 
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Fig. 6. Photometric (Ve) residuals on the first night of the caustic crossing, relative 
to the linear limb-darkening fit, with a parabolic ri{t) shown in Fig. |^. The individual 
residuals are extremely small, but the trend is highly significant, showing that the real 
light curve peaks earlier than the model by about 0.1 days. 
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Fig. 7. Difference in {Ve — Ie) color between the two highly-magnified caustic images 
and the source as a whole, as a function of r}, the separation of the source center from 
the caustic in units of the source radius. The color shift is extremely small and only 
marginally significant on the first two nights, but rises to 0.046 ±0.001 on the third night 
and 0.27 ± 0.02 on the last, when the caustic images are equivalent to looking at the 
outer 5.5% of the source. Since the source is a K3 giant, this limb has the color of a MO 
giant. 



